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Abstract 

l  - 

The  steady-state  dynamic  fields  near  a  rapidly  propagating  crack  tip 
in  an  elastic  perfectly-plastic  material  have  been  investigated  for  the  case 
of  Mode-Ill  fracture.  For  arbitrary  values  of  the  dimensionless  crack-tip 
speed  (M)  the  inner  solution  consists  of  a  central-fan  field  ahead  of  the 
crack  tip  and  a  uniform  field  in  its  wake.  It  is  shown  that  the  inner  solution 
is  valid  in  a  "boundary  layer"  which  shrinks  on  the  crack  tip  in  the  limit  of 
vanishing  M.  For  small  M  the  outer  solution  was  found  as  a  regular  perturbation 
expansion  in  M,  with  the  quasi-static  solution  as  its  first  term.  A  uniform 
expansion  over  the  polar  angle  8  measured  from  the  plane  of  the  crack  was  also 
obtained;  its  first  term  displays  the  connection  between  the  inner  and  outer 
solutions . 
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1.  Introduction 

In  a  recent  paper  the  authors  have  investigated  dynamic  effects  on  the 
fields  of  stress  and  strain  near  a  rapidly  propagating  crack  tip  in  an  elastic 
perfectly-plastic  material  [1].  For  the  steady-state  case  it  was  found  that 
the  dynamic  near-tip  fields  can  be  expressed  as  simple-wave  solutions  of  the 
governing  system  of  hyperbolic  partial  differential  equations.  These  solu¬ 
tions  are  independent  of  the  dimensionless  distance  to  the  crack  tip,  r/r^  , 
but  they  do  depend  specifically  on  the  dimensionless  crack-tip  speed  M.  For 
Mode-Ill  crack  propagation  the  simple-wave  solution  in  the  near-tip  field  is 
a  combination  of  a  centered-fan  field  and  a  uniform  field.  Explicit 
expressions  have  been  presented  in  Refs.[l]  and  [2], 

The  solutions  that  were  obtained  in  Ref.[l]  show  some  anomalies  in  the 
transition  from  the  dynamic  to  the  quasi-static  solution.  As  the  crack-tip 
speed,  M,  decreases  the  expressions  for  the  stresses  reduce  to  the  ones  for 
the  corresponding  quasi-static  solution,  as  might  be  expected  on  the  basis  of 
intuitive  reasoning.  This  is  however  not  true  for  the  strains,  which  become 
unbounded  in  the  limit  of  vanishing  M.  In  Ref.[l]  it  was  speculated  that  the 
transition  from  dynamic  to  quasi-static  conditions  with  decreasing  crack-tip 
speed  is  effected  because  the  dynamic  solution  is  asymptotically  valid  in  a 
small  edge  zone,  which  shrinks  on  the  crack  tip  in  the  limit  of  vanishing 
crack-tip  speed. 

In  this  paper  this  non-uniform  transition  has  been  investigated  in  detail 
in  the  plane  of  the  crack  for  the  case  of  crack  propagation  in  anti-plane 
strain  (Mode-Ill) .  It  is  shown  that  for  small  crack-tip  speeds  the  complete 
near-tip  solution  consists  of  the  outer  solution,  which  is  a  regular 
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perturbation  expansion  in  M,  with  the  quasi-static  solution  as  the  first  term, 
and  the  inner  solution,  which  is  of  completely  different  nature  with  a  strong 
influence  of  dynamic  effects.  The  inner  solution  is  valid  in  an  edge-zone  which 
is  analogous  to  a  boundary  layer.  The  first  term  in  a  uniform  expansion  over  the 
polar  angle  0  measured  from  the  plane  of  the  crack  displays  the  connection  between 
the  inner  and  outer  solutions. 

The  governing  equations  and  the  boundary  conditions  are  stated  in  Section  2. 
The  outer  solution,  its  relation  to  the  quasi-static  solution,  as  well  as  its 
inadequacy  in  the  immediate  vicinity  of  the  crack  tip  have  been  discussed  in 
Section  3.  The  results  for  the  inner  solution  have  been  summarized  in  Section 
4.  It  is  shown  in  Section  5,  that  certain  functions  which  define  the  fields  for 
small  polar  angle  from  the  plane  of  the  crack,  satisfy  a  system  of  coupled 
ordinary  differential  equations.  An  investigation  of  the  singular  points  in 
the  phase  plane  and  an  inspection  of  the  trajectories  of  the  solutions  of  these 
equations  as  the  distance  to  the  crack  tip  decreases,  reveals  that  the  expressions 
of  Section  4  do  indeed  provide  the  solution  in  an  edge  zone.  Exact  solutions  in 
implicit  form  for  the  coupled  ordinary  differential  equations  mentioned  above 
have  also  been  obtained  in  Section  5,  and  these  solutions  reproduce  the  inner 
solution  as  well  as  the  outer  solution.  Finally  the  coupled  equations  have 
been  solved  numerically,  and  the  results  have  been  plotted  vs.  r/r^  and  M,  to 
show  the  transition  of  the  inner  solution  in  the  edge  zone  to  the  outer  solution 
in  the  "far-field". 
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2.  Governing  Equations 

A  coordinate  system  (x,y,z)  is  attached  to  the  moving  crack  tip  as  shown 
in  Fig.  1.  The  assumption  that  a  steady-state  has  been  established  relative  to 
the  moving  crack  tip,  implies  that  absolute  time  derivatives  take  the  following 
forms  relative  to  the  moving  coordinate  system: 


V  > 


-v  ax(  ); 


3tt<  > 


v23  (  ) 


XX 


In  the  moving  coordinate  system  steady-state  deformation  in  anti-plane  strain 
is  defined  by  a  displacement  in  the  z-direction,  which  is  denoted  by  w(x,y) . 
The  corresponding  stress  components  are  a^^x.y)  and  a^^(x, y) . 

In  view  of  (2.1b)  the  equation  of  motion  becomes 


3  a  +  8  a 
x  xz  y  yz 


pv23  w 

XX 
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(2.1a,b) 


The  Tresca  yield  condition  is 

a2  +  a2  *  k2  '  , 
xz  yz 

where  k  is  the  yield  stress  in  shear.  By  virtue  of  (2.1a),  the  Prandtl-Reuss 
equations  for  an  elastic  perfectly-plastic  solid  reduce  to  the  forms 

A.>  0 

Here u  is  the  elastic  shear  modulus,  and  A  is  a  non-negative  proportionality 
factor,  which  may  vary  in  space. 

In  the  region  of  plastic  deformation  the  yield  condition  is  identically 
satisfied  by 

a  ■  -  k  sinw  ,  a  ■  k  cosw 

xz  yz 

Elimination  of  A  from  (2.4)  and  (2.5),  and  the  use  of  (2.2)  and  (2.6a,b)  then 
yields  the  following  system  of  equations 


(2.2) 


(2.3) 


(2.4) 

(2.5) 


(2.6a,b) 
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cosu  3  «  +  sinai  3  w  +  M2  3  oi  »  0 
x  y  k  x  x 


costo  3  w  +  sinoi  3w  +  —  3  in  ■  0 
xx  y  x  vi  x 


where 

w  ■  3  w  ,  and  w  =  3  w 
xx  y  y 

M  -  v/(u/p)^ 

The  condition  of  vanishing  o^z  on  the  crack  faces  yields  by  virtue  of  (2.6b) 
y  **  0,  x  <  0:  a)  *  ir/2 

Displacement  anti- symmetry  relative  to  y  ■  0  implies  that  w  ■  0  for  y  *  0. 
Hence  *  0,  and  thus  by  (2.6a) 
y  *  0,  x  >  0:  w  »  0 

In  the  polar  coordinate  system  shown  in  Fig.  1,  (2.7)  and  (2.8)  have 
the  forms 


(2.7) 

(2.8) 

(2.9a, 

(2.10) 

(2.11) 


(2.12) 


cos(ui-0)  3  u  +  —  sin(u>-9)  3-ui  + 
r  r  d 


+  (cos9  3  w  -  sin0  3„w )  -  0 
k  r  x  r  9  x 


(2.13) 


cos(id-0)  3  w  +  —  sin(w-0)  3„w  + 
r  x  r  v  x 


,  k  (cos03  oi - sin0  3aoi)  ■  0 

+  —  r  r  y 

v 


3 .  Far-Field  Solution 


Let  us  consider  regular  series  expansions  for  u  and  with  respect  to 

the  "Mach  number"  M  defined  by  (2.10) 

a)  ■  oi  +  Mm,  +  .  .  . 
o  1 

w  ■  w  +Mw,  +  ... 

X  xo  xl 


(2.14) 


(3.1) 


(3.2) 
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For  the  leading  terms  n>  and  w^o  we  obtain  the  following  equations 

cos(w  -0)  3  (d  +  sin(o)  -0)3ato  ■  0 
o  r  o  r  o  oo 


(3.3) 


cos(u)  -0)  3  w  +  —  sin (a  -0)  3aw  + 
o  r  xo  r  o  0  xo 


k.  1 

+  —  (cos0  3  i») - sin0  3au)  ) 

U  r  o  r  0o 


m  *  9 
o 


(3.4) 


These  equations  govern  the  quasi-static  problem,  which  was  discussed  by 
Chitaley  and  McClintock  [  3  ].  For  small  values  of  r/rp  (where  rp  defines  the 
boundary  of  the  plastic  domain)  we  obtain  (see  Ref.  3  ) 


(3.5) 


wxo  *  "u  sin0  8,11  + 

P 

where  f(0)  is  an  arbitrary  function  of  9. 

It  is  of  interest  to  examine  the  magnitude  of  the  Inertia  term,  M23xwx  » 
which  would  correspond  to  the  quasi-static  solution.  By  the  use  of  (3.5)  and 
(3.6)  we  find 

M23  w  =  jr—  — 2  [l-in(—  )  -  ^21]  sin29 
XX  2p  r  rp  cos  0J 


(3.6) 


(3.7) 


On  the  other  hand,  if  the  stress  derivatives  appearing  in  Eq.(2.2)  are  computed 
on  the  basis  of  (3.5)  and  (3.6),  it  is  found  that  they  vanish  identically. 
Equation  (3.7)  then  suggests  that  the  quasi-static  approximation  is  not  valid 
in  a  small  neighborhood  of  the  crack-tip  defined  by  r/rp  ~  0(M2)  (we  shall 
make  a  more  accurate  estimate  in  a  later  Section).  The  subsequent  terms  in  the 
series  (3.1)  and  (3.2)  do  not  remove  this  nonuniformity  of  the  approximation. 
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Thus  the  regular  expansion  with  respect  to  M  cannot  be  accepted  in  the 
immediate  vicinity  of  the  crack-tip.  It  appears  that  this  solution  represents 
an  outer  or  "far-field"  expansion. 

4,  Near-Field  Solution 

In  a  recent  paper  Achenbach  and  Bunayevsky  [  1  ]  have  shown  that  the 
solution  (oi,wx)  to  Eqs.(2.7)  and  (2.8)  is  either  singular  near  the  crack  tip, 
or  it  is  represented  by  a  centered  fan-field  in  combination  with  a  uniform 
field.  Since  it  was  not  possible  to  obtain  a  singular  solution,  the  centered 
fan  +  uniform  field  solution  was  considered  in  some  detail.  The  solution  was 
obtained  in  the  following  form  f  1  ,  see  also  2  ] : 

0  1  0  1  e* 

w  cos  ^  [Msin20  +  (l~M2sin0)^cos0]  (4.1) 

x 

a  *  -k[ (i-M2sin20)^  -  Mcos0]  sin0  (4.2) 

xz 

a  =  k[ (1-M2sin20) *  cos0  +  Msin20]  (4.3) 

yz 

* 

0  <  0  <  n 

a  *  -k  ,  a  *0,  w  ■  -kir/pM  '  (4.4a, 

xz  yz  x 

where 

9  *  tan  ^(1/M)  (4.5) 

In  the  limit  M  0  the  stress  fields  reduce  to  the  corresponding  quasi¬ 
static  fields.  This  is,  however,  not  true  for  the  strains,  which  become 
unbounded  as  M  -*•  0.  Difficulties  of  this  kind  in  the  transition  from  the 
dynamic  to  the  quasi-static  solution  as  M  +  0  could  have  been  anticipated 
from  the  structure  of  the  governing  equations.  As  M  +  0  the  two  distinct 
families  of  characteristic  curves  of  Eqs.(2.6)  and  (2.7),  which  are  defined 


by 


dx  M 

—  =  COSO)  ±  — ; - 

dy  sitro 

merge  into  one  family,  defined  by 
dx 


dy 


COS(0 


Such  a  degeneracy  usually  leads  to  a  non-uniform  transition,  and  the 
appearance  of  a  "boundary  layer".  Some  examples  are  given  by  Cole  [4  ]. 

Another  indication  of  a  non-uniform  transition  is  given  by  the  form  of 
the  equations  in  the  hodograph  plane.  To  transform  (2.7)  and  (2.8)  to  the 
hodograph  plane  we  introduce  the  following  changes  of  variables  [  5J. 

.  .  3w  .  .  3w 

3x  _  ,-l  _ x  _  _  j-1  _ x 

9(0  9y  *  9(0  9x 


3x  _  ,-1  dm  dy  _  ,- 

3W  ~  J  >  a...  _  J 


3y 


9w 


1  1“ 
9x 


where  J  is  the  Jacobian 


3  9w  3w  . 

j  _  _9(o  _ x _ x  dm 

3x  9y  9x  9y 

In  the  hodograph  plane  Eqs.(2.7)  and  (2.8)  take  the  form 


9y  .  9x  „?  k  9y 

cosio  t* - sinoo  - - M4  — =  0 

3w  9w  p  9(0 

X  X 


-cosor^-  +  sinio  -r—  +  —  =  0 

9(0  3to  p  3w 

x 


This  linear  system  of  equations  can  be  reduced  to  a  single  equation  of  the 
second  order 

(V  if*  .  „2  +  M2cotio  -  -±-  |*-  -  0 

P  »  9  ^  9  3(o  sinui  3w 


3w2 

x 


3(0  2 


(4.6) 

(4.7) 

(4.8a,h) 

(4.9a,b) 

(4.10) 

(4.11) 

(4.12) 


(4.13) 
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When  M  -*•  0  the  second  order  derivative  32y/3w2  disappears,  which  indicates 
that  an  asymptotic  expansion  with  respect  to  M  cannot  be  uniform. 

On  the  basis  of  the  foregoing  observations  it  may  be  assumed  that  the 
field  near  the  crack  tip  is  of  different  forms  in  two  zones.  The  outer 
solution  in  the  "far  field"  is  a  regular  expansion  in  M  as  given  by  (3.1)  and 
(3.2),  where  the  leading  term  represents  the  quasi-static  solution.  The  effect 
of  inertia  is  relatively  small  in  this  far-field.  Inertia  is,  however, 
important  in  the  inner  solution  in  the  near  field,  represented  by  (4.1)-(4.5). 
The  near-field  can  be  thought  of  as  a  "boundary  layer".  In  the  present 
geometry  the  terminology  "edge  zone"  is,  however,  more  appropriate.  In  the 
edge-zone  the  inertia  effects  appear  to  remove  the  singularity  of  the  quasi¬ 
static  strain  component. 

In  the  next  section  we  consider  the  matching  of  the  inner  solution  and  the 
outer  solution,  and  the  shrinking  of  the  edge-zone  on  the  crack  tip  as  M  +  0. 


5.  Solution  in  the  Plane  of  the  Crack 

Matching  of  the  far-field  and  near-field  solutions  has  not  been  accomplished 
for  arbitrary  values  of  the  polar  angle  0.  It  is,  however,  possible  to  investi¬ 
gate  the  matching  for  small  values  of  0,  by  seeking  a  solution  in  the  following 
form 

u)  =  a)^(r)  0  +  oi^(r)  0^  +  .  .  .  (5.1) 

w^  =  w_^(r)  ®  +  wx3^r^  ®3  +  •  •  •  (5.2) 

Substitution  of  (5.1)  into  (2.3)  and  (2.4),  and  collecting  the  terms  of  equal 
powers  of  0  ,  gives  for  the  first  approximation  the  following  equations 


+ 


^(w.j-1) 

r 


dw 


xl 


dr 


(5.3) 


dwxl  ,  wxl  .  k 

dr  r  v  'dr 


-)  =  0 


At  this  stage  it  is  useful  to  introduce  new  variables 


J  =  Mw  +  —  a)  , 
+  xl  y  1 


„  k 

Mw  . - u), 

xl  p  1 


Inserting  (5.5)  into  (5.3)  and  (5.4)  we  arrive  at 


dJ. 


««,)  am  ^  <vj-W+  • 0 


(1-M)  ~  +  (1-M)J_-^  (J+-J_)J_  =  0 

where  the  following  change  of  variables  has  been  used 
a  -  -Hn(r/rp) 

Let  us  consider  the  solutions  to  Eqs.(5.6)  and  (5.7)  in  the  phase  plane 


(5.4) 


(5.5ar 


(5.6) 

(5.7) 


(5.8) 
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(J+,J  ).  The  singular  points  are  the  solutions  to  the  following  equations 

(14*1) J+  -  ~  (J+-J_)J+  =  °  (5.9) 

( 1-M) J_  -  ^  =  0  (5.10) 

Equations  (5.9)  and  (5.10)  define  three  singular  points,  whose  position  and 
character  are  given  by 


1) 

J+  =  J- 

=  0  ; 

stable  focus 

(5.11)  j 

2) 

J_  =  0; 

J+  =  (14*0  &  ; 

unstable  node 

(5.12) 

3) 

<-* 

+ 

II 

o 

j_  =  - ( 1-M)  ^  ; 

saddle  point 

(5.13) 

The  nature  of  the  singular  points  given  by  (5. 11) -(5. 13)  determines  the  phase 
flow  pattern,  i.e.,  the  trajectories  of  the  solutions  as  a  varies,  shown  in 
Fig.  2  (see  e.g.  Ref.  [6]). 

Since  no  limit  cycle  or  center  point  exists,  the  trajectories  beneath  the 
separatrix  line  AA'  run  towards  infinity  as  a  -*■  «  (r+O) .  Along  these  lines 
J+  and  J  ,  and  consequently  and  w^,  become  unbounded.  However  infinite 
growth  of  leads  to  an  oscillation  of  the  stress  field  (with  increasing 
frequency)  which  is  not  admissible  from  the  physical  point  of  view.  Thus  this 
domain  of  the  phase  plane  falls  outside  our  consideration. 

The  phase  flow  above  the  separatrix  line  AA'  tends  to  the  origin  J+  *  J_  *  0 
so  that  for  a  ->•  °°  (r+0)  we  have  J+  -*■  0,  J_  -*  0,  and  consequently  -*•  0, 
w^  -*■  0.  This  result  would  correspond  to  a  uniform  zero  field  ahead  of  the 
crack  tip.  The  hyperbolicity  of  the  governing  equations  (2.7)  and  (2.8)  then 
would  imply  a  zero  field  for  the  whole  loading  domain,  which  does,  however,  not 


J 
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satisfy  the  boundary  condition  on  the  crack  faces  and  therefore  has  to  be 
rejected. 


Thus  we  have  no  other  alternative  but  the  trajectory  along  the 
separatrix  AA'  .  There  is  however  a  "forbidden"  region  in  the  phase  plane 
which  is  defined  by  the  condition  A  ^  0.  It  follows  from  (2.4),  (2.5)  and 
(5.1),  (5.2)  that 

w  i 

k  A  *  -  — 


The  requirement  A  _>  0  leads  to  the  inequality 


xl 


2M  <J++J-)  ±  0 


( 


which  implies  that  we  have  to  consider  the  part  of  the  separatrix  which  lies 
to  the  left  of  the  bisector  J+  =  -J  . 

Thus,  when  a  -*•  °°  (r+0)  the  solution  moves  along  the  line  AA'  toward  the 
point  3.  Let  us  now  find  the  solution  in  the  neighborhood  of  the  point  3. 
Linearization  of  Eqs.(5.6)  and  (5.7)  near  this  point  leads  to  the  following 
equations 

dJ+  ^  _  2M_ 
da  =  1+M  J+ 


~  -  -  J ,  +  J  +  (1-M)  — 
da  +  -  u 

In  terms  of  the  variable  r  the  solution  has  the  form 
-  c  if!  <r/rp>2"/<1'm> 

J  ■  -  (1-M)  —  +  C(r/r  )2M^1+M) 

up 

where  C  is  an  arbitrary  constant  and  (5.8)  has  been  used.  Hence  by  virtue  of 
(5.5)  we  obtain 


(5.14) 

(5.15) 


(5.16) 

(5.17) 

(5.18) 

(5.19) 


12 


»  .  -  -  \  (l-M)  +  C  i2»  (r/r 
xl  M  4M  p 

Hence  we  have 


(5.20) 

(5.21) 


ui  s  w  0  -  (l-M)e  +  0  [ (r/r 

1  P 


(5.22) 


w  ;  w  0  =*  -  (l-M)e  +  0  [(r/r  )2M/(1+M>  } 
x  xl  yM  p 

The  dominant  terms  of  the  expressions  (5.22)  and  (5.23)  coincide  with  those  of 
the  expansion  of  the  near  field,  (4.1)  and  (4.4),  for  small  0.  This  result 
confirms  the  validity  of  the  simple  wave  solution  near  the  crack  tip. 

Equations  (5.6)  and  (5.7)  can  be  integrated  to  yield  the  following 


(5.23) 


solution  in  implicit  form 

J  -  *  Cr/r  )2M/<1W,) 
+  P 


lJJ 


(l-M)  /  (1+M) 


(5.24) 


F(1^,  1  +  1i,  -B<f> 

p 


2M/(1+M) 


-2M/(1-M) 


-1 


)-  "  l^M  2k  J-  +  B  J-(f  )  (5'25) 


Here  A  and  B  are  constants  of  integrations  and  F(p,q,r,s)  is  the  hypergeometrical 
function.  The  trajectory  along  the  separatrix  AA'  corresponds  to  B  ■  0.  This 
follows  from  the  observation  that  the  separatrix  crosses  the  line  J+  ■  0  at  the 
point  3,  which  is  defined  by  (5.13) 

It  can  now  be  verified  that  for  (r/r^)  -*■  0,  (5.24)  and  (5.25)  yield 
J+  -*■  0  and  J_  -*•  (2k/u)(l-M),  which  correspond  to  ui^  ■*>  l-M  and  wx^  -*•  -(k/yM)  (l-M) . 
These  results  agree  with  (5.20)  and  (5.21). 

Next  we  consider  (5.24)  and  (5.25)  in  the  limit  M  -*■  0.  It  is  noted  from 
(5.5)  that  for  M  ■  0  we  have  J+  -  -J  .  Since  J  <0,  the  constant  A  in  (5.24) 
must  equal  unity:  A  ■  1.  For  M  ^  0,  the  combination  of  (5.5a,b)  and  (5.24) 
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then  yields 


*n(Mwxi  +  v u) "  Hi  *n  (f->  +  m *n  |Mwxi "  £  “ 


l-M 


To  expand  (5.26)  for  small  M  we  require  r/r  >  5  where  5  is  such  that 

P 


-  (f->  -  0(|) 


Equation  (5.26)  then  yields 

k.  .  f  /T  \  .  /k 


wxl  =  -jju  [fcn(^  )  -  bi  in(-^u)]  +  0(M) 


In  deriving  this  expression,  it  was  taken  into  account  that  J  <  0. 

For  the  hypergeometrical  function  in  Eq.(5.25)  we  find  for  small  M 


F(1-  1  +  ii,  -  <f> 

P 


2M/1+M  -2M/1+M 


lJ_l 


) 


1  +  Sr  ~  i)  +  -  g5  e#*>  +  °(m2> 


2M  ,  *  '  2M  , 
Here  B(x)  is  the  g-f unction 


o  /  \  a  _1  «  _ n^_^ _  1 

X  2  n=tQ  x(x+l) . . .  (x+k)  2n 


By  virtue  of  (5.29),  (5.30)  and  (5.25)  we  obtain  for  M-K) 


4  +  0(M)  -  -  (1+M)  (Mw  .  -  -  oa) 

2  2k  xl  y 


For  (r/r^)  >  5  the  formulas  (5.26)  and  (5.31)  imply 


U)  -  1  ,  w  -  -  to  (— )  +  -  ftnC— ) 

1  ’  xl  u  r  v  u 

P 

This  result  agrees  with  the  quasi-static  solution  given  by  (3.5)-(3.6),  which 
coincides  with  the  dominant  terms  of  the  outer  solution  for  small  9,  as 
given  by  (3.1)— (3.2). 


(5.26) 


(5.27)  ; 


(5.28) 


(5.29) 


(5.30) 


(5.31) 


(5.32) 


1 


To  follow  in  some  detail  the  phase  flow  when  M  +  0  it  is  convenient  to 


consider  the  (w^co)  phase  plane  shown  in  Fig.  3.  The  point  of  intersection  of 
the  separatrix  AA'  with  the  co-axis  should  be  located  between  1-M  and  14M.  The 
angle  y  of  intersection  of  the  J  -  and  J  -lines,  is  given  by 


y  =  ff  -  2  tan  ^(1/M), 

which  tends  to  zero  when  M  0.  The  abscissa  of  the  saddle  point  3  is 

=  1-M  and  the  abscissa  of  the  node  2  is  u>^  =  1+M,  which  both  tend  to  1  as 
M  ■>  0.  Thus  the  singular  points  2  and  3  shift  to  infinity  as  M  -*•  0.  For  M  =  0 
the  phase  flow  has  degenerated  as  shown  in  Fig.  4.  The  separatrix  has  turned 
into  a  straight  line  which  is  parallel  to  the  w^-axis  and  has  the  abscissa 
ui-1.  The  magnitude  w^  has  become  unbounded  along  this  line. 

In  summary,  for  small  values  of  the  angle  0  the  solution  given  by  (5.24) 
and  (5.25)  matches  the  near  field  (4.1)  and  (4.3)  to  the  far  field  (3.5)  -  (3.6) 
for  small  values  of  the  angle  0  .  The  region  defined  by  in(r/r^)  ~  0(M  ^)  can 
be  identified  as  an  edge  zone. 

The  implicit  form  of  (5.24)  and  (5.25)  Is  not  convenient  for  the  computation 
of  curves.  Equations  (5.6)  and  (5.7)  have,  therefore,  been  solved  numerically. 
The  function 


vxl  -  (J+  +  JJ/2M 

has  been  plotted  in  Fig.  5.  The  corresponding  quasi-static  solution,  which  is 
indicated  by  M  -  0,  has  also  been  plotted.  The  quasi-static  solution  for 
is  singular  at  r/r^  *  0,  while  the  dynamic  solution  remains  bounded.  The 
curve  for  M  -  0.01  is  very  close  to  the  one  for  M  =  0,  i.e.,  to  in 

the  region  r/r^  >  6,  where  6  is  the  "length"  of  the  edge  zone  in  the  plane  of 
the  crack.  An  estimate  of  the  length  of  the  edge  zone  is  indicated  in  Fig.  5. 
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